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Abstract — Low-density parity check (LDPC) codes are a signif- 
icant class of classical codes with many applications. Several good 
LDPC codes have been constructed using random, algebraic, and 
finite geometries approaches, with containing cycles of length at 
least six in their Tanner graphs. However, it is impossible to 
design a self-orthogonal parity check matrix of an LDPC code 
without introducing cycles of length four. 

In this paper, a new class of quantum LDPC codes based 
on lines and points of finite geometries is constructed. The 
parity check matrices of these codes are adapted to be self- 
orthogonal with containing only one cycle of length four in 
each pair of two rows. Also, the column and row weights, and 
bounds on the minimum distance of these codes are given. As 
a consequence, these codes can be encoded using shift-register 
encoding algorithms and can be decoded using iterative decoding 
algorithms over various quantum depolarizing channels. 

I. Introduction 

Low density parity check (LDPC) codes are a capacity- 
approaching (Shannon limit) class of codes that were first 
described in a seminal work by Gallager [8]. In Tanner [21], 
LDPC codes were rediscovered and presented in a graphical 
interpretation (codes over graphs). Iterative decoding of LDPC 
and turbo codes highlighted the importance of these classes of 
codes for communication and storage channels. Furthermore, 
they have been used extensively in many applications [7], [15], 
[16]. 

There have been several notable attempts to construct regu- 
lar and irregular good LDPC codes using algebraic combina- 
torics and random constructions, see [16], [20], and references 
therein. Liva et al. [16] presented a survey of the previous 
work done on algebraic constructions of LDPC codes based 
on finite geometries, elements of finite fields, and RS codes. 
Furthermore, a good construction of LDPC codes should have 
a girth of the Tanner graph, of at least six [15], [16]. 

Quantum information is sensitive to noise and needs error 
correction, control, and recovery strategies. Quantum block 
and convolutional codes are means to protect quantum infor- 
mation against noise and decoherence. A well-known class of 
quantum codes is called stabilize codes, in which it can be 
easily constructed using self-orthogonal (or dual-containing) 
classical codes, see [3], [5], [11] and references therein. 
Recently, subsystem codes combine the features of decoher- 
ence free subspaces, noiseless subsystems, and quantum error- 
correcting codes, see [2], [4], [13], [14] and references therein. 

Quantum block LDPC codes have been proposed in [17], 
[18]. MacKay et al. in [17] constructed sparse graph quantum 



LDPC codes based on cyclic matrices and using a computer 
search. Recently, Camera el al. derived quantum LDPC codes 
in an analytical method [6]. Hagiwara and Imai constructed 
quasi-cyclic (QC) LDPC codes and derived a family of 
quantum QC LDPC codes from a nested pair of classical 
codes [10]. 

In this paper, we construct LDPC codes based on finite 
geometry. We show that the constructed LDPC codes have 
quasi-cyclic structure and their parity check matrices can 
be adapted to satisfy the self-orthogonal (or dual-containing) 
conditions. The motivations for this work are that (i) LDPC 
codes constructed from finite geometries can be encoded using 
linear shift-registers. The column weights remain fixed with 
the increase in number of rows and length of the code, (ii) The 
adapted parity check matrix has exactly one cycle with length 
four between any two rows and many cycles with length of at 
least six. (iii) A class of quantum LDPC codes is constructed 
that can be decoded using known iterative decoding algorithms 
over quantum depolarizing channels; some of these algorithms 
are stated in [19]. 

Notation: Let q be a prime power p and F q be a finite field 
with q elements. Any two binary vectors v = (vi, V2, ■ ■ ■ , v n ) 
and u = (iti, 1*2, . . . ,u n ) are orthogonal if their inner product 
vanishes, i.e., 2~^™=i ViUi m °d 2 = 0. Let H be a parity check 
matrix defined over F2, then H is self-orthogonal if the inner 
product between any two arbitrary rows of H vanishes. 

II. LDPC Code Constructions and Finite 
Geometries 

A. LDPC Codes 

Definition 1: An (p, A) regular LDPC code is defined by a 
sparse binary parity check matrix H satisfying the following 
properties. 

i) p is the number of one's in a column. 

ii) A is the number of one's in a row. 

iii) Any two rows have at most one nonzero element in 
common. The code does not have cycles of length four 
in its Tanner graph. 

iv) p and A are small in comparison to the number of rows 
and length of the code. In addition, rows of the matrix H 
are not necessarily linearly independent. 

The third condition guarantees that iterative decoding algo- 
rithms such as sum-product or message passing perform well 
over communication channels. In general it is hard to design 



regular LDPC satisfying the above conditions, see [15], [16], 
[20] and references therein. 

B. Finite Geometry 

Finite geometries can be classified into Euclidean and 
projective geometry over finite fields. Finite geometries codes 
are an important class of cyclic and quasi-cyclic codes because 
their encoder algorithms can be implemented using linear 
feedback shift registers and their decoder algorithms can 
be implemented using various decoding algorithms such as 
majority logic (MLG), sum-product (SPA), and weighted BF, 
see [12], [15], [16]. 

Definition 2: A finite geometry with a set 
{p 1 ,p 2 ,...,p n }, a set of I lines {L 1 ,L 2 ,... 
integer pair (A, p) is defined as follows: 

i) Every line Li passes through p points. 

ii) Every point pi lies in A lines, i.e., every point p 
intersected by A lines. 

iii) Any two points p\ and pj can define one and only one 
line Lk in between. 

iv) Any two lines Li and Lj either intersect at only one point 
Pi or they are parallel. 

Therefore, we can form a binary matrix H = [hi j] of size 
I x n over F2. The rows and columns of H correspond the / 
lines and n points in the Euclidean geometry, respectively. If 
the if/2 line Li passes through the point p. 
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Fig. 1. Euclidean geometry with points n = 4 and lines ( = 6 

geometry with n = 4, I = 6, A = 3, and p = 2. We 
can construct the incidence matrix H based on this geometry 
where every point and line correspond to a column and row, 
respectively. For p « I and A << n, The matrix H is a 
sparse low density parity check matrix. In this example, the 
matrix Heg-i is given by 
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We call the Euclidean geometry defined in this type as a 
Type-I EG. The Tanner graph of Type-I EG is a regular 
bipartite graph with 11 code variable vertices and I check-sum 





(a) 



Fig. 2. (a) Euclidean geometry with n = 4 points (check-sum) and I = 6 
lines (code variables) (b) This Tanner graph corresponds to a self-orthogonal 
parity check matrix. 



vertices. Also, each variable bit vertex has degree A and each 
check-sum has degree p. 

If we can take the transpose of this matrix Heg-i, then we 
can also define a (p, A) LDPC code with length I and minimum 
distance is at least p + 1. The codes defined in this type are 
called LDPC codes based on Type-II EG. In this type, any 
two rows intersect at exactly one position. 

C. Adapting the Matrix Heg-ii to be Self-orthogonal 

Let Heg-ii be a parity check matrix of a regular LDPC 
code constructed based on Type-II EG Euclidean geometry. 
We can construct a self-orthogonal matrix H^'L/j from 
Heg-ii m two cases. 

Case 1. If the number of one's in a row is odd and any two 
rows intersect at exactly one position, i.e., any line connects 
two points. As shown in Fig. 2, the Tanner graph corresponds 
to a self-orthogonal parity check matrix H°E G h _ n if and only 
if every check-sum has even degree and any any two check- 
sum nodes meet at even code variable nodes. This condition 
is the same as every row in the parity check matrix ^°eg-ii 
has an even weight and any two rows overlap in even nonzero 
positions. 



H 



EG-II 



(2) 



The vector 1 of length n is added as the last column in 



TTOrth 
"-EG-II- 

Case 2. Assume the number of one's in a line is even 
and any two rows intersect at exactly one position. We can 
construct a self-orthogonal parity check matrix I^eg-ii as 
follows. We add the vector 1 along with the identity matrix I 
of size n x n. We guarantee that any two rows of the matrix 
H°EG h -n intersect at two nonzero positions and every row has 
an even weight. 
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D. Characteristic Vectors and Matrices 
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Let n be a positive integer such that n = q m — 1, where 
m = ord„(g) is the multiplicative order of q modulo n. Let a 
denote a fixed primitive element of F g ™ . Define a map z from 
F* m to such that all entries of z(a l ) are equal to except 
at position i, where it is equal to 1. For example, z(a 2 ) = 
(0, 1,0,..., 0). We call z(a k ) the location (or characteristic) 
vector of or. We can define the location vector z(a J+:,+1 ) 
as the right cyclic shift of the location vector z(a* +J ), for 
< j < n — 1, and the power is taken module n. The location 
vector can be extended to two or more nonzero positions, for 
example, the location vector of a 2 , a 3 and a 5 is given by 
z(a 2 ,a 3 ,a 5 ) = (0, 1, 1, 0, 1, 0, . . . , 0). 

Definition 3: We can define a map A that associates to an 
element F* m a circulant matrix in F 2 l x " by 



(q m -l)/(q-l), 



(6) 



A(a*) = 



z(a l+1 ) 



\ z(a l+n - r ) ) 



(4) 



By construction, A(a k ) contains a 1 in every row and column. 

We will use the map A to associate to a parity check matrix 
H = (hij) in (F* m ) the (larger and binary) parity check 
matrix H = (A(hij)) in F£ xn . The matrices A(hij)'s are 
n x n circulant permutation matrices based on some primitive 
elements hij as shown in Definition 3. 

III. Constructing Self-Orthogonal Cyclic LDPC 
Codes from Euclidean Geometry 

In this section we construct self-orthogonal algebraic Low 
Density Parity Check (LDPC) codes based on finite geome- 
tries. Particulary, there are two important classes of finite 
geometries: Euclidean and projective geometry. 

A. Euclidean Geometry EG(rn, q) 

We construct regular LDPC codes based on lines and points 
of Euclidean geometry. The class we derive has a cyclic 
structure, so it is called cyclic LDPC codes. Cyclic LDPC 
codes can be defined by a sparse parity check matrix or by a 
generator polynomial and can be encoded using shift-register. 
Furthermore, they can be decoded using well-known iterative 
decoding algorithms [15], [16]. 

Let q be power of a prime p, i.e. q = p s for some 
integer s > 2. Let EG(m, q) be the m-dimensional Euclidean 
geometry over F g for some integer m > 2. It consists of 
p-ms _ qm p 0m ]- s an( j ever y point is represented by an m- 
tuple, see [12]. A line in EG(m,q) can be described by a 
1 -dimensional subspace of the vector space of all m-tuples 
over F q or a coset of it. The number of lines in EG(m, q) is 
given by 



(g— l)(^-l)/(g-l), 



(5) 



and each line passes through q points. Every line has q( m x > — 
1 lines parallel to it. Also, for any point in EG(m, q), there 



lines intersect at this point. Two lines can intersect at only one 
point or they are parallel. 

Let F q m be the extension field of F q . We can represent each 
element in F g ™ as an m-tuple over F g . Every element in the 
finite field F q ™ can be looked as a point in the Euclidean 
geometry EG{m, q), henceforth F q m can be regarded as the 
Euclidean geometry EG(m, q). 

Let a be a primitive element of F q ™. q m points of 
EG(m, q) can be represented by elements of the set 
{0, 1, a, a 2 , . . . , a q ~ 2 }. We can also define a line L as the 
set of points of the form {a + 7 b | 7 S F q }, where a and b 
are linearly independent over F q . For a given point a, there 
are (q m — l)/(q— 1) lines in EG(m, q) that intersect at a. 

Type-I EG. Let 11 = q m — 1 be the number of points 
excluding the original point in EG(m,q). Assume L be 
a line not passing through 0. We can define the binary vector 



v L = (vi,vi, . . .,v n ), 



(7) 



where Vi = 1 if the point a 1 lies in a line L. The vector 
is called the incidence vector of L. Elements of the vector 
Vi correspond to the elements 1, a, a 2 , ... , a™ -1 . aL is also 
a line in EG(m, q), therefore is a right cyclic-shift of 
the vector vl. Clearly, the lines L, aL, . . . , a n ~ 1 L are all 
different. But, they may not be linearly independent. 

Consider the vectors Li,aLi, . . . ,a n ~ 1 Li. We can con- 
struct an n x n matrix H; in the form 



Hi = 



( ^ 

a\ Li 



\ 



(8) 



Clearly, Hi is a circulant matrix with column and row weights 
equals to q, the number of points that lie in a line a? Li, for 
< j < n — 1. Hi has size of n x n. The total number of 
lines in EG(m, q) that do not pass through the origin are 
given by 

(9) 

They can be partitioned into (q" 1 ^ 1 — l)/(g— 1) cyclic classes, 
see [16]. Every class Hi can be defined by an incidence vector 



Li as {Li, aL u a 2 Li, ...,oi 
l)/(q-l). Let 1 < £ < Of"" 1 
defined as 



n_1 ii} for 1 < i < (<f 



I)/ (q - 1), then Heg,i is 



T~Leg,. 



Hi H2 



H 



(10) 



For each cyclic class Hi, we can form the matrix H; over 
F2 of size n x n. Therefore, is a circulant binary matrix 
of row and column weights of q. 

If we assume that there are 1 < I < (g m_1 — l)/(q — 1) 
incidence lines in EG(m, q) not passing through the origin, 
then we can form the binary matrix 



H 



EG, I 



Hi H 2 



H, 



(11) 



The matrix Heg,j£ consists of a £ sub-matrices H; of size 
nxn and it has column and row weights £q and q, respectively. 
The null space of the matrix Heg,£ gives a cyclic EG-LDPC 
code of length n = q rn — 1 and minimum distance £q + 1, 
whose Tanner graph has a girth of at least six, see [16], [20]. 

The Tanner graph of Type-I EG is a regular bipartite graph 
with q m — 1 code variable vertices and I check-sum vertices. 
Also, Each variable bit vertex has degree p = q and each 
check-sum has degree A = £ q. 

Type-II EG. We can take the transpose of the parity check 
matrix 1~t(EG,e) over Fg m as defined in Type-I to define a new 
parity check matrix with the following properties, see [12]. 



H 



EG,. 



= H{ 



Hi 



(12) 



So, the matrix Hf is the transpose matrix of Hi- Consequently, 
we can define the binary matrix Heg/ 



a EGj 



H 



(13) 



Let £ = (g m_1 - l)/(g - 1), then the matrix H EG e has the 
following properties 

i) The total number of columns is given by £n = (q™ 1 ^ 1 — 
l)(g--l)/(g-l). 

ii) Number of rows is given by n — q m — 1. 

iii) The rows of this matrix correspond to the nonorigin points 
of EG(m, q) and the columns correspond to the lines in 
EG(m, q) that do not pass through the origin. 

iv) A = £q = qiq™- 1 -l)/{q-l) = (q m - l)/(q- 1) - 1 is 
the row weight for £ = (q m_1 — l)/(g — 1). Also p = q 
is the column weight. 

v) Any two rows in H|J G e have exactly one nonzero element 
in common. Also, any two columns have at most one 
nonzero element in common. 

vi) The binary sub-matrix Hj has size (q™ — 1) x (q m — 1). 
Also, it can be constructed using only one vector vl that 
will be cyclically shifted q m — 1 times. 

B. QC LDPC Codes 

The matrix H|J G t defines a quasi-cyclic (QC) LDPC code 
of length N = £n'= (gC™" 1 ) - l)(q m - l)/(q - 1) for £ = 



(q" 



!)/(« 



1). The matrix H 



EG,. 



has n 



1 



iii) Every two columns intersect at most at one nonzero 
position. Every two rows have exactly one and only one 
nonzero position in common. 
From this definition, the minimum distance of the LDPC 
code defined by the null-space of H|J G e is at least p+1. This 
is because we can add at least p+1 columns in the parity check 
matrix H E g i to obtain the zero column (rank of H E g i i s a * 
least (p + 1)). Furthermore, the girth of the Tanner graph for 
this matrix H, is at least six, see [7], [20]. This is a (p, A) QC 
LDPC code based on Type-II EG. 



C. Self-orthogonal QC LDPC Codes 

We can define a self-orthogonal parity check matrix H^ G ' l £ 
from Type-II EG construction as follows. The binary matrix 
h egj of size nx £nfov I < £ < (g™' 1 - l)/(q - 1) has 
row and column weights of A = £q and p = q, respectively. 
Let 1 be the column vector of size (q m — 1) x 1 defined as 
1 = (1,1,..., 1) T . If the weight of a row in H|J G e is odd, 



then we can add the vector 1 to form the matrix H 



rth 
EG,i 



H 



T 

EG.l 



1 . Also, if the weight of a row in H BG e is even, 
then we can add the vector 1 along with the identity matrix 



orth 



of size (q m - 1) x (q m - 1) to form H^ G ' 



H 



EG, 



I 



Therefore, we can prove that VL°eq\ is self-orthogonal as 



shown in the following Lemma. 

Lemma 4: The parity check matrix H^ G l £ defined as 



H 



orth 
EG J 



bJ I l 



, for odd 



H 



ill , for even £q 



is self-orthogonal. 

Proof: From the construction Type-II EG, any two dif- 
ferent rows intersect (overlap) in exactly one nonzero position. 
If £q is odd, then adding the column vector 1 will result an 
even overlap as well as rows of even weights. Therefore, the 
inner product mod 2 of any arbitrary rows vanishes. Also, 
if £q is even, adding the columns 1 | I will produce row of 
even weights and the inner product mod 2 of any arbitrary 
rows vanishes. ■ 



H 



orth 
EG,£ 



rows that are not necessarily independent. We can define a QC 

LDPC code over F 2 as the null-space of the matrix H| G t of the P 3 ^ check matrix H sG> is at least 1 + 1 



has size n x N for odd £q where n = q m — 1, 
N = n£ + 1, and 1 < £ < (q^- 1 ) - l)/(q - 1). Also, it has 
length N = n(l + 1) + 1 for even £q. 

The minimum distance of the LDPC codes constructed in 
this type can be shown using the BCH bound as stated in the 
following result. 

Lemma 5: The minimum distance of an LDPC defined by 

orth 



sparse circulant sub-matrices of equal size. The matrix H EG e 
with parameters (p, A) has the following properties. 

i) p = q is the weight of a column c,;. p does not depend 
on m, hence length of the code can be increased without 
increasing the column weight. 

ii) A = Iq is the weight of a row r^. A depends on m, but 
the length of the code increases much faster than A. 



IV. Quantum LDPC Block Codes 

In this section we derive a family of LDPC stabilizer codes 
derived from LDPC codes based on finite geometries. Let P = 
{I, X, Z,Y = iXZ} be a set of Pauli matrices defined as 



/ = 



1 ) > X ~ ( 1 ) > Z ~ ( -1 



(14) 



and the matrix Y is the combination of the matrices X bit 

/ —i 

flip and Z phase-flip defined as Y = iXZ = I 
Clearly, 

= z 2 = r 2 = l 



A well-known method to construct quantum codes is by 
using the stabilizer formalism, see for example [1], [5], [9], 
[17] and references therein. Assume we have a stabilizer group 
S generated by a set {Si, S2, ■ ■ ■ , <Sn-fc} such that every two 
row operators commute with each other. The error operator Sj 
is a tensor product of n Pauli matrices. 



Sj = E 1 ®E 2 



E n , Ei G P. 



Sj can be seen as a binary vector of length 2n [5], [17]. 
A quantum code Q is defined as +1 joint eigenstates of the 
stabilizer S. Therefore, a codeword state \ip) belongs to the 
code Q if and only if 



Sj\ip) = \ip) for all Sj e S. 



(15) 



CSS Construction: Let G and H be two binary matrices 
define the classical code C and dual code C , respectively. 
The CSS construction assumes that the stabilizer subgroup 
(matrix) can be written as 



H 





(16) 



where H and G are k x n matrixes satisfying HG = 0. 
The quantum code with stabilizer S is able to encode n — 2k 
logical qubits into n physical qubits. If G = H, then the self- 
orthgonality or dual-containing condition becomes HH T = 0. 
If C is a code that has a parity check matrix H, then C x C C. 

Constructing Dual-containing LDPC Codes: Let us con- 
struct the stabilizer matrix 



S 



stab 



The matrix H%ae 

is a binary self-orthogonal matrix as 
shown in Section III-C. We replace every nonzero element 
in H'eq'i by the Pauli matrix X to form the matrix Hx- 
Similarly, we replace every nonzero element in H ^q\ by the 
Pauli matrix Z to form the matrix Hz- Therefore the matrix 
S s tab is also self-orthogonal. We can assume that the matrix 
Hx corrects the bit-flip errors, while the matrix Hz corrects 
the phase-flip errors, see [1], [17]. 

Lemma 6: A quantum LDPC code Q with rate (n — 2k)/n 
is a code whose stabilizer matrix S s tab of size 2k x 2n has a 
pair (p, A) where p is the number of non-zero error operators 
in a column and A is the number of non-zero error operators 
in a row. Furthermore, S s t a b is constructed from a binary self- 
orthogonal parity check matrix H^.'^ of size k x n. 

Using Lemma 6 and LDPC codes given by the parity check 
matrix as shown in Section III-C, we can derive a class 

of quantum LDPC codes as stated in the following Lemma. 

Theorem 7: Let H^/^ be a parity check matrix of an 
LDPC code based on EG{m, q), where n = q m — 1 and 











(17) 



1 < £ < (q m 1 — l)/(q — 1). Then, there exists a quantum 
LDPC code Q with parameters \[N, N — 2n, > q + 1]]2 where 
N = £n+ 1 for odd £q and N = (£+ l)n + 1 for even £q. 

Proof: By Lemma 4, H^Tq 1 ^ is self-orthogonal. Using 
Lemma 6, there exists a quantum LDPC code with the given 
parameters. ■ 

V. Conclusion 

We constructed a class of quantum LDPC codes derived 
from finite geometries. The constructed codes have high rates 
and their minimum distances are bounded. They only have one 
cycle of length four between any two rows and many cycles 
of length of at least six. A new class of quantum LDPC codes 
based on projective geometries can be driven in a similar way. 
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